 &MODEL
 NLIST   =           1,
 NSTEP   =        2000,
 EPSI    =  1.0000000E-07,
 WGHT    =   1.000000    ,
 NOBS    =        1000,
 NVAXZ   =           4,
 NVAX    =           3,
 NORDER  =           0,
 NVAZ    =           1,
 NLAMZ   =           1,
 NLAMY   =           1,
 NLAMX   =           3,
 NSCALE  =           1,
 NEXPY   =           0,
 NRES    =           0,
 NCOR    =           0,
 NCOV    =           0,
 VARNAME = BC4HG    ,
 NVARIANT        =          17
 /
 &HETERO
 NUMZ    =           4, 9*0,
 NUMLZ   =           1, 9*0,
 DELTA   = 10*0.000000000000000E+000  ,
 NTDEL   =           1, 9*0,
 VLZ     = 10*0.000000000000000E+000  ,
 NTLZ    =           1, 9*0
 /
 &YBOXCOX
 VLY     =  0.000000000000000E+000,
 NTLY    =           1,
 NUMLXY  =           0
 /
 &XBOXCOX
 NUMX    =           1,           2,           3, 47*0,
 NUMLX   =           1,           2,           3, 47*0,
 VLX     = 10*0.000000000000000E+000  ,
 NTLX    = 3*1, 7*0
 /
 &YXZDATA
 YXZFILE = TijL2.in                                                                        ,
 YXZFMT  = (5e15.7)                                                                        
 /

VARIABLE       MINIMUM         MEAN          MAXIMUM   NULL OBS.      SCALE
DEPENDENT
Tij         0.5899500D+04  0.2573965D+08  0.3419374D+09     0     0.1000000D-07
INDEPENDENT
Ti          0.1580900D+05  0.1860662D+06  0.3513510D+06     0     0.1000000D-04
Tj          0.6962000D+04  0.1851223D+06  0.4337820D+06     0     0.1000000D-04
Util_BC4    0.3868082D-07  0.7346229D-03  0.1399558D-01     0     0.1000000D+02
CONSTANT    0.1000000D+01  0.1000000D+01  0.1000000D+01     0     0.1000000D+01
HETEROSKEDASTICITY
TiTj        0.6827108D+09  0.3439253D+11  0.1093478D+12     0     0.1000000D-10


MRS of X(l) with respect to X(k) - at the sample means
===================================================================================================================================
        X(k) Ti          Tj          Util_BC4    CONSTANT  
X(l)
===================================================================================================================================

Ti           0.1000D+01  0.1229D+01  0.2776D+09  0.3413D+06

Tj           0.8137D+00  0.1000D+01  0.2259D+09  0.2777D+06

Util_BC4     0.3602D-08  0.4426D-08  0.1000D+01  0.1229D-02

CONSTANT     0.2930D-05  0.3601D-05  0.8135D+03  0.1000D+01


INITIAL ESTIMATES

             LAM    PARAMETER      STD-ERROR        STUDENT-T    GRADIENT       DERIV.Y         ELAST.Y
                                                                                DERIV.E(Y)      ELAST.E(Y)


Ti            X1  0.5325132D+00  0.6598434D-01      8.0702964  0.2799840D-11  0.7366572D+02      0.5325132
                                                                              0.1169046D+00      0.5325132

Tj            X2  0.6510872D+00  0.5906770D-01     11.0227274  0.1252758D-11  0.9052803D+02      0.6510872
                                                                              0.1475469D+00      0.6510872

Util_BC4      X3  0.5837007D+00  0.1616134D-01     36.1171024 -0.3518445D-10  0.2045165D+11      0.5837007
                                                                              0.1703872D+03      0.5837007

CONSTANT      X0  0.9767329D+00  0.1192295D+00      8.1920401  0.7217138D-11  0.2514076D+08      0.9767329
                                                                              0.3300557D+00      0.9767329

VAR(W)            0.1263045D+01  0.5648508D-01     22.3606798  0.0000000D+00

INITIAL VALUES
 
TiTj       Z1  0.0000000D+00
LAM(Z1)        0.0000000D+00
 
LAM(Y)         0.0000000D+00
LAM(X1)        0.0000000D+00
LAM(X2)        0.0000000D+00
LAM(X3)        0.0000000D+00


LOG-LIKELIHOOD =  0.9344055D+03

NUMBER OF OBSERVATIONS              =           1000

PEARSON-R2(E) (unadjusted)          =         0.2435
PEARSON-R2(E) (adjusted for D.F.)   =         0.2412

PSEUDO-R2(L)  (unadjusted)          =         0.5780
PSEUDO-R2(L)  (adjusted for D.F.)   =         0.5767

PROBABILITY of Y to be at the limit =         0.0000
MEAN VALUE of E(Y)                  =  0.3379181D+00

MEAN of residual W                  =  0.9464949D-14
VARIANCE of residual W              =  0.1263045D+01
95% CONFIDENCE INTERVAL for residual W
Lower Bound = -0.2202751D+01   Upper Bound =  0.2202751D+01




initialization cards for control routine follow - -  id card - no. parameters , no. steps , step size 
minimization tolerance - starting values of parameters - relative uncertainties of parameters (error estimates) -  
and starting direction  (relative parameter increments for 1st step)
identification line
         6      2000   0.1E+01   0.1E-06
  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00
  0.10000E+01  0.10000E+01  0.10000E+01  0.10000E+01  0.10000E+01  0.10000E+01
  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00

first entry to fcn , parameters of current step follow - - f,x(i)
  -0.9344054864D+03  0.00000D+00  0.00000D+00  0.00000D+00  0.00000D+00  0.00000D+00  0.00000D+00

new absolute minimum , parameters for fcn entry     5  follow - - f,x(i)
  -0.9453577071D+03 -0.20584D-01  0.00000D+00  0.63533D-01 -0.11152D-02  0.55658D-02 -0.20584D-01

new absolute minimum , parameters for fcn entry     8  follow - - f,x(i)
  -0.9511145538D+03 -0.12198D+00 -0.28518D-02  0.83242D-01 -0.66550D-02  0.34882D-01  0.23987D-01

new absolute minimum , parameters for fcn entry    11  follow - - f,x(i)
  -0.9619195649D+03 -0.36112D+00 -0.17323D-01  0.58541D-01 -0.20913D-01  0.11509D+00 -0.17235D-01

new absolute minimum , parameters for fcn entry    15  follow - - f,x(i)
  -0.9773609044D+03 -0.44444D+00  0.37645D+00  0.42075D-01 -0.17528D+00  0.12556D+01 -0.41653D-01

new absolute minimum , parameters for fcn entry    18  follow - - f,x(i)
  -0.9800622875D+03 -0.40781D+00  0.13364D+00  0.52654D-01 -0.23962D+00  0.15217D+01 -0.26841D-01

new absolute minimum , parameters for fcn entry    21  follow - - f,x(i)
  -0.9802932421D+03 -0.41533D+00  0.15035D+00  0.53684D-01 -0.42481D+00  0.15178D+01 -0.27884D-01

new absolute minimum , parameters for fcn entry    24  follow - - f,x(i)
  -0.9817284034D+03 -0.55165D+00  0.32274D+00  0.73880D-01 -0.32093D+00  0.17030D+01 -0.34811D-01

new absolute minimum , parameters for fcn entry    28  follow - - f,x(i)
  -0.9835982486D+03 -0.89101D+00  0.67672D+00  0.67449D-01 -0.32765D+00  0.18465D+01 -0.21315D-01

new absolute minimum , parameters for fcn entry    31  follow - - f,x(i)
  -0.9852243966D+03 -0.12647D+01  0.80764D+00  0.70009D-01 -0.32383D+00  0.16020D+01 -0.34328D-01

new absolute minimum , parameters for fcn entry    34  follow - - f,x(i)
  -0.9858228364D+03 -0.16414D+01  0.11340D+01  0.72440D-01 -0.27907D+00  0.16056D+01 -0.29355D-01

new absolute minimum , parameters for fcn entry    37  follow - - f,x(i)
  -0.9862208565D+03 -0.19865D+01  0.12549D+01  0.72497D-01 -0.33804D+00  0.16356D+01 -0.22261D-01

new absolute minimum , parameters for fcn entry    40  follow - - f,x(i)
  -0.9865878685D+03 -0.25545D+01  0.15121D+01  0.75061D-01 -0.31056D+00  0.16135D+01 -0.31792D-01

new absolute minimum , parameters for fcn entry    43  follow - - f,x(i)
  -0.9867338665D+03 -0.30797D+01  0.16903D+01  0.74309D-01 -0.32042D+00  0.16579D+01 -0.29651D-01

new absolute minimum , parameters for fcn entry    45  follow - - f,x(i)
  -0.9867671222D+03 -0.33886D+01  0.18203D+01  0.75865D-01 -0.32869D+00  0.16298D+01 -0.31362D-01

new absolute minimum , parameters for fcn entry    47  follow - - f,x(i)
  -0.9867731824D+03 -0.34747D+01  0.18348D+01  0.76192D-01 -0.33085D+00  0.16262D+01 -0.31767D-01

new absolute minimum , parameters for fcn entry    49  follow - - f,x(i)
  -0.9867741005D+03 -0.35389D+01  0.18593D+01  0.76372D-01 -0.32975D+00  0.16314D+01 -0.31688D-01

new absolute minimum , parameters for fcn entry    51  follow - - f,x(i)
  -0.9867741171D+03 -0.35443D+01  0.18610D+01  0.76317D-01 -0.33002D+00  0.16308D+01 -0.31711D-01

new absolute minimum , parameters for fcn entry    52  follow - - f,x(i)
  -0.9867741173D+03 -0.35447D+01  0.18611D+01  0.76317D-01 -0.33009D+00  0.16309D+01 -0.31718D-01

new absolute minimum , parameters for fcn entry    53  follow - - f,x(i)
  -0.9867741173D+03 -0.35447D+01  0.18611D+01  0.76317D-01 -0.33008D+00  0.16309D+01 -0.31717D-01

new absolute minimum , parameters for fcn entry    54  follow - - f,x(i)
  -0.9867741173D+03 -0.35447D+01  0.18611D+01  0.76317D-01 -0.33008D+00  0.16309D+01 -0.31717D-01


the function (fcn) has not improved by    0.10000D-06 after 3 steps which include a minimum - - -


therefore the run has been terminated
  -0.9867741173D+03 -0.35447D+01  0.18611D+01  0.76317D-01 -0.33008D+00  0.16309D+01 -0.31717D-01


all numbers of last entry refer to final pt. of run


the lowest value of the function (fcn) was  -0.9867741173D+03 for entry   54 - - x(i) follows
 -0.35447E+01  0.18611E+01  0.76317E-01 -0.33008E+00  0.16309E+01 -0.31717E-01



run statistics follow - -

 nfcn=   58 nst=   21


MRS of X(l) with respect to X(k) - at the sample means
===================================================================================================================================
        X(k) Ti          Tj          Util_BC4    CONSTANT  
X(l)
===================================================================================================================================

Ti           0.1000D+01  0.2037D+01  0.3233D+09 -0.2454D+08

Tj           0.4910D+00  0.1000D+01  0.1587D+09 -0.1205D+08

Util_BC4     0.3093D-08  0.6299D-08  0.1000D+01 -0.7590D-01

CONSTANT    -0.4075D-07 -0.8299D-07 -0.1317D+02  0.1000D+01


FINAL ESTIMATES

             LAM    PARAMETER      STD-ERROR        STUDENT-T      STUDENT-T    GRADIENT       DERIV.Y         ELAST.Y
                                   (CONDIT.)        (CONDIT.)      (CONDIT.)                   DERIV.E(Y)      ELAST.E(Y)
                                                    (PAR = 0)      (PAR = 1)

Ti            X1  0.8024562D+02  0.3267204D+03      0.2456094                -0.1250843D-06  0.5500485D+02      0.3976179
                                 0.1037223D+02      7.7365844                                0.8149687D+02      0.4234748

Tj            X2  0.7606232D-08  0.2643593D-07      0.2877233                -0.1475476D+02  0.1120305D+03      0.8057354
                                 0.5985419D-09     12.7079349                                0.1109985D+03      0.9045117

Util_BC4      X3  0.1484748D+01  0.5232567D+00      2.8375137                 0.4659576D-06  0.1778419D+11      0.5075700
                                 0.3635835D-01     40.8365116                                0.1744056D+12      0.5495196

CONSTANT      X0 -0.1928592D+03  0.7351382D+03     -0.2623442                -0.4204714D-07 -0.1349852D+10    -52.4425088
                                 0.3074795D+02     -6.2722633                               -0.1344468D+10    -52.8562628

VAR(W)            0.1828117D+02  0.1045432D+02      1.7486704                 0.0000000D+00
                                 0.1093625D+01     16.7161161
 
TiTj      (E) Z1 -0.1195985D-19  0.1446205D-18     -0.0826981                 0.1519379D+15  0.2550222D-04      0.0340753
                                 0.1710732D-20     -6.9910709                               -0.1340034D-03     -0.2178829
LAM(Z1)   (E)     0.1861077D+01  0.4917315D+00      3.7847424      1.7511123 -0.4450336D-04
 
LAM(Y)    (E)     0.7631728D-01  0.1766434D-01      4.3204137    -52.2908011  0.1227749D-04
LAM(X1)   (E)    -0.3300774D+00  0.3534529D+00     -0.9338655     -3.7630968 -0.2731702D-04
LAM(X2)   (E)     0.1630877D+01  0.2826975D+00      5.7689838      2.2316342 -0.2259683D-05
LAM(X3)   (E)    -0.3171737D-01  0.1348871D-01     -2.3514011    -76.4874625  0.1918217D-05


LOG-LIKELIHOOD = -0.1743391D+05

NUMBER OF OBSERVATIONS              =           1000

PEARSON-R2(E) (unadjusted)          =         0.3125
PEARSON-R2(E) (adjusted for D.F.)   =         0.3063

PSEUDO-R2(L)  (unadjusted)          =         0.8844
PSEUDO-R2(L)  (adjusted for D.F.)   =         0.8834

PROBABILITY of Y to be at the limit =         0.0000
MEAN VALUE of E(Y)                  =  0.2750202D+08

MEAN of residual W                  = -0.2230045D-01
VARIANCE of residual W              =  0.1828117D+02
95% CONFIDENCE INTERVAL for residual W
Lower Bound = -0.8402571D+01   Upper Bound =  0.8357970D+01


CPU TIME:      0.39 sec
